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We show that the analytic iV-soliton solution of the Camassa-Holm (CH) shallow-water 
model equation converges to the nonanalytic iV-peakon solution of the dispersionless CH 
equation when the dispersion parameter tends to zero. To demonstrate this, we develop a 
novel limiting procedure and apply it to the parametric representation for the iV-soliton 
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as well as various identities among the Hankel determinants to facilitate the asymptotic 
analysis. We also provide a new representation of the iV-peakon solution in terms of the 
Hankel determinants. 

KEYWORDS: Camassa-Holm equation, soliton, peakon, parametric representation 



* E-mail address: matsuno@yamaguchi-u.ac.jp 



1 



1. Introduction 

The Camassa-Holm (CH) equation is a model equation describing the unidirectional 
propagation of nonlinear shallow-water waves over a flat bottom. 1 " 4 ) It may be written 
in a dimensionless form as 

u t + 2k 2 u x - u txx + 3uu x = 2u x u xx + uu xxx , (1.1) 

where u = u(x,t) is the fluid velocity, k is a positive parameter related to the phase 
velocity of the linear dispersive wave and the subscripts t and x appended to u denote 
partial differentiation. Like the Korteweg-de Vries equation in shallow-water theory, the 
CH equation is a completely integrable nonlinear partial differential equation with a rich 
mathematical structure. Because of this fact, a large number of works have been devoted 
to the study of the equation from both physical and mathematical points of view. We do 
not provide a comprehensive review on various consequences established for the equation. 
Instead, we refer to an excellent survey given in ref. 5). Our main concern here is 
solutions to the CH equation for both k ^ and k — 0. Although equation (1.1) can 
be transformed to the corresponding equation with k = by a Galilean transformation 
x' = x + K 2 t, t' — t,u' — u + k 2 , solutions under the same vanishing boundary condition 
at infinity have quite different characteristics. In fact, in the case of k ^ 0, it typically 
exhibits analytic multisoliton solutions (the so-called iV-soliton solution with N being an 
arbitrary positive integer) like the usual soliton equations. 6 ~ 15 ) We also remark that it 
has singular cusp soliton solutions as well as solutions consisting of an arbitrary number 
of solitons and cusp solitons. 16 " 19 ) When k, — 0, on the other hand, eq. (1.1) becomes 
a dispersionless version of the original CH equation. It admits a new kind of solitary 
waves which have a discontinuous slope at their crest. 1 ' 2 -* For this unique feature, they 
are now termed peakons. Hence, the peakon has a nonanalytic nature unlike the smooth 
soliton. One of its remarkable properties is that the motion of peakons is described by 
a finite dimensional completely integrable dynamical system. This fact was used in refs. 
20 and 21 to construct the general iV-peakon solution which represents the interaction 
of iV peakons. The careful inspection of the interaction of peakons reveals that it occurs 
elastically in pairs and the total effect of the collision is the phase shift whose explicit 
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formula has been given in a closed form. 2 ' 9,21 ) See also a recent paper concerning the 
detailed investigation of the dynamics of two peakons. 22 ) An important issue about the 
iV-peakon solution is how one can reduce it from the analytic iV-soliton solution by taking 
the singular limit k — > 0. Although this problem has received considerable attention, it 
has been resolved only partially. Indeed, in the case of N = 1, the convergence of the single 
soliton to the single peakon has been demonstrated using the explicit form of the 1-soliton 
solution. 2 ' 5 ^ Also, an analysis without recourse to the explicit form of the soliton solution 
has been carried out in a general context by employing an abstract theory of dynamical 
systems. 23 ' 24 ^ However, the treatment of the general iV-soliton case has remained open. 
Quite recently, the case N = 2 was completed by two different ways. One will be given in 
ref. 25 while the other is included here as a special version of the general iV-soliton case. 

The purpose of this paper is to demonstrate the convergence of the iV-soliton solution 
to the iV-peakon solution for the general N by using the explicit iV-soliton solution of the 
CH equation. In §2, we present a parametric representation for the iV-soliton solution 
of the CH equation which offers a relevant form to the subsequent asymptotic analysis. 
In §3, the convergence of the 1-soliton solution to the 1-peakon solution is exemplified to 
explain the new idea used in the limiting procedure. In §4, we perform the corresponding 
limit for the iV-soliton solution and show that it reproduces the iV-peakon solution given 
by Beals et al. 21 ) We also obtain a new representation of the 7V-peakon solution in terms 
of the Hankel determinants. Section 5 is devoted to the concluding remarks. In Appendix 
A, we give a proof of the formula which enables us to rewrite the tau-functions for the N- 
soliton solution of the CH equation in terms of the Hankel determinants. In Appendix B, 
we establish various identities among the Hankel determinants which are used effectively 
to simplify the limiting waveform of the 7V-soliton solution. 

2. The iV-Soliton Solution of the CH Equation 

The iV-soliton solution of the CH equation (1.1) may be represented in a parametric 



form 11 ' 15 ) 




(2.1a) 



x(y,t) = ^ + ln ^ + d, 



(2.16) 



3 



with 



fi= ex P 

(t=0,l 

= exp 

n=o,i 



N 



,i=l l<i<j<N 
' N 



i=l 



l<i<j<N 



where 



ii = h(y- KCit - y i0 ) , c, 



2k 2 



(i = l,2,...,N), 



g r. _ 



1 + ftfcj ' 



e 7ij = 



l-/t 2 A; 2 ' 
(0</cfci<l), (i = l,2,...,iV), 

= l,2,...,JV;i^ j). 



(2.2a) 
(2.26) 

(2.3a) 
(2.36) 
(2.3c) 



(hi + /cj) 2 

Here, fcj and y i0 are soliton parameters characterizing the amplitude and phase of the ith 
soliton, respectively, q is the velocity of the ith soliton in the (x, t) coordinate system and d 
is an integration constant. We assume that q > and q ^ Cj for i ^ j = 1, 2, N). 
In (2.2), the notation ^ 01 implies the summation over all possible combination of 
/ii = 0, l,/i2 = 0, 1, ...,/In = 0, 1. The following coordinate transformation (x,t) — > (y,t') 
has been introduced to parametrize the iV-soliton solution 



dy = rdx — urdt, dt' = dt, 



(2.4a) 



where the variable r is defined by the relation 



2 , 2 

r = u — u xx + k . 



(2.46) 



Note in (2.1) that the time variable t' is identified with the original time variable t by 
virtue of (2.4a). 

The 1-soliton solution is given by 

2K 2 c x k\ 



u(y,t) 



1 + K 2 kj + (1 - K 2 k 2 ) COsh^ 



x(y,t) = - + In 

K 



(1 - Kki)e^ + 1 + nki 
((1 + Kki)e^ + 1 - Kki 



+ d. 



(2.5a) 
(2.56) 
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It represents a solitary wave travelling to the right with the amplitude c\ — 2k 2 and velocity 
c\. The property of the 1-soliton solutions has been explored in detail as well as its peakon 
limit. See, for example ref. 5. For completeness, we write the explicit 1-peakon solution 
of equation (1.1) with k — 0. It reads 

u(x,t) = Cl e-\ x - Clt - Xl °\, (2.6) 

where x w represents the initial position of the peakon. In view of the invariance property 
of the dispersionless version of the CH equation under the transformation x' — x, t' — 
—t, v! = —u, it also admits a peakon solution (2.6) with a negative amplitude. It represents 
a peakon with depression propagating to the left. 

3. The Peakon Limit of the 1-Soliton Solution 

We first consider the limiting procedure for the 1-soliton solution. This will be helpful 
to understand the basic idea in developing the procedure for the general iV-soliton solution. 
We start with the tau-functions for the 1-soliton solution which are the most important 
constituents in our analysis. These are given by (2.2) with N = 1. Before proceeding to 
the limit, we find it convenient to shift the phase constant y±o as yio — > y±o + 4>i/k\ or in 
terms of the phase variable £i — > £i — <pi so that 

h = l + e«\ (3.1a) 

/ 2 = l + ^ e fc (3.16) 

where we have put v\ = e _<?il . The limit k — > is taken in such a way that the amplitudes 
of the soliton and peakon given respectively by C\ — 2k 2 and C\ coincide and remain 
finite. 5 ^ This can be attained by taking the limit k — * with the soliton velocity c\ in the 
(x,t) coordinate system being fixed. We see from (2.3a) that the appropriate limit of the 
wavenumber k\ is carried out by taking nk\ — > 1. The limiting procedure described here 
may be called the peakon limit. It now follows from (2.3) that various wave parameters 
have the following leading-order asymptotics in the peakon limit 

K 2 K 2 Ai 

nki ~ 1 , ui ~ — = — k 2 , (3.2a) 

ci 2ci 4 v 1 
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e ?i _ e Kfci(^-cit-xio) ^ h e x- Cl t-x 10 ^ (3.26) 

where Ai = 2/ci and rr 10 = Vio/k. In passing to the last line of (3.2b), we have used (2.1b) 
to eliminate the y variable and the constant d has been absorbed in the phase constant 
x w . Substituting (3.2b) into (3.1), the leading terms of f\ and / 2 are found to be as 

fi ~ 1 + T^i, (3.3a) 
h 

/ 2 ~l + e 2 A 2 ^, (3.36) 

J2 

where we have put 

zi = e x - Clt - Xl °, e = ^i ( 3 - 3c ) 
for simplicity. If we introduce the new quantity / by / = /2//1, we can deduce from (3.3) 
that 

7 "TmT- (3 ' 4) 

The expression (3.4) yields the following quadratic equation for /: 

/ 2 + (^-l)/-e 2 A 2 z 1 + O( e 3 ) = 0. (3.5) 

If we solve this equation, we can express / as a function of x and t. At this instant, it 
is crucial to observe that / > since both f\ and f 2 are positive quantities. A positive 
solution is then substituted into (2.1a) to obtain u. In performing the differentiation, 
however, we must replace the t derivative d/dt by d/dt + ud/dx in accordance with the 
coordinate transformation (2.4a). With this notice in mind, we can rewrite (2.1a) as 

«=(» + «|)ln / . (3.6) 

Solving (3.6) with respect to u to express it in terms of a single variable / and its deriva- 
tives, we have 

u = jk-. (3.T) 

The derivatives ft and f x in (3.7) are obtained simply if one differentiates (3.5) by t and 
x, respectively. To be more specific, as e — > 

ft ~ 27 + ^-1 ' (3 - 8a) 
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- 2f + Zl -l ' ^ 
where we have used the relations zi )t = — ci^i and zi jX = z\ which are derived from (3.3c). 

After a manipulation using (3.5) and (3.8), we arrive at the expression of u in terms of /: 

u ~ - 7T — . (3.9) 

The final step of the limiting process is to solve (3.5) under the condition / > and 
then take the limit k — > after substituting a positive solution into (3.9). Although the 
analytical expression is obtained for / by quadrature, we need only the series solution. 
This fact is crucial in developing the peakon limit of the general iV-soliton solution where 
equation corresponding to (3.5) becomes an algebraic equation of degree N + 1 whose 
analytical solution is in general not available. 
Now, we expand / in powers of e as 

f = /( 0) + e/ (l) + e 2 /( 2) + (310) 

and insert this expression into (3.5). Comparing the coefficients of e n (n = 0,1,...), we 
obtain a system of algebraic equations for f^ n \ the first two of which read 

/(°) 2 + ( Zl - = 0, (3.11a) 

2/(o) / (D + {zi _ !)/(!) = o, (3.116) 
_ _ X \ Zl = 0. (3.11c) 

The above equations can be solved immediately to obtain positive solutions.. Indeed, if 
Zi < 1 ( x — C]t — x w < 0), then 

/ ~ /(°) = 1 - ^. (3.12a) 
Substitution of this result into (3.9) yields 

u ~ cizi = Cl e x - Clt - Xl °. (3.126) 
If, on the other hand, Z\ > l(or x — c±t — x w > 0), then 

/ (0) = =0, / ~ /(V = \\z l /{z l - l)e 2 , (3.13a) 



u ~ ci^f 1 = cie- (x - cl * _a!lo) . (3.136) 

It follows from (3.12) and (3.13) that in the limit k — > (or equivalently e — > by (3.3c)) 
w has a limiting waveform 

u = Cie -l*-cit-*iol. (3.14) 
This expression is just the 1-peakon solution (2.6) of the CH equation with k — 0. 

4. The Peakon Limit of the iV-Soliton Solution 

4-1 Formulas for determinants 

The peakon limit of the general iV-soliton solution can be taken along the lines of 
the 1-soliton case. However, the calculation involved is quite formidable. To perform the 
calculation in an effective manner, we first define the following determinants which are 
closely related to the iV-peakon solution: 
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i ••• i 


in) = 




^2 " - \ n 




\ n—1 


a"- 1 ••• Ar 1 




n 





l<l<m<n 



A r , 



m+l 



A 



m+l 



A 



m+n— 1 



Here 



A. 

A 

■"-m+n 
N 



A 



m+n— l 
A 

^m+n 



A 



2(n-l)+m 



(4.1) 



(4.2a) 



Ei 



1=1 

K Ci 



1,2, ...,7V), 



(4.26) 



(4.2c) 



where in (4.2) m is an arbitrary integer and n is a nonnegative integer less than or equal to 
N. For n greater than TV, M m) = 0. We use the convention A^) = 1, = 1, L>J m) = 
A m . The quantity A„ is the square of the Vandermonde determinant whereas D n m ^ is the 
determinant of a symmetric matrix. It is a Hankel determinant. We use some properties 
of Hankelians in the following analysis. 
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Let DT'(ii,i2, ip,ji, 32, -,j q ) (h < h < ••• < i P , ji < h < ••• < j q , 1 < P,q < N) 
be a determinant which is obtained from by deleting rows ii,i 2 , ■■■■,i p and columns 
ji, j 2 , ■■■,jq, respectively. Then, the following Jacobi formula holds which will play a 
central role in the present analysis: 26 ) 

D%\D%\(l,n + 2;l,n + 2)=D^ 

(4.3) 

By virtue of the definition (4.2), we see that 

D™(l,n + 2;l,n + 2) = D$»+ 2 \ (4.4a) 

D^ 2 (l;l) = Dt + r\ (4-46) 

D%\(n + 2;n + 2) = D%\, (4.4c) 

£>£> (1; n + 2) = £>£> (n + 2; 1) = (4.4d) 
Hence, (4.3) can be rewritten in the form 

nM n( m + 2 ) — n( m+2 )n( m ) f n( m +i)\ 2 (a k\ 

LJ n+2 LJ n — U n+\ U n+\ ~ [^n+l J ■ 

The determinant has an alternative expression in the form of a finite sum 21 ' 

D M = J- A n (tut2,...,in)(K\ l2 ..X n rE tl E i2 ...E ln , (n = 1,2, ...,N). (4.6) 

i<ii<i 2 <...<i„<iV 

This formula is very useful in rewriting the tau- functions. Note that (n = 1, 2, N) 
are positive definite since A n > and > (i — 1, 2, iV). We give a simple proof of 
(4.6) in Appendix A. 

4-2 The peakon limit of the N-soliton solution 

We first shift the phase variables as £j — > £j — = 1,2, ...,N) in (2.2) and take 
the peakon limit /t/cj — > 1 with fixed Cj(i = 1,2, ...,N). Using (3.2) and the asymptotic 
e 7lJ ~ (Aj — Aj) 2 /t 4 /16, the leading-order asymptotics of the tau-functions fx and f 2 can 
be written in the form 

fi ~ 1 + ( ^ &n(ii,i2,..,i n )z h z i2 ...z in , (4.7a) 

n=l V/2/ i<i 1<i2< ... <in <N 
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f 2 „ i + e n(n+1) If) (Ai 1 A ia ...Ai B ) 2 A B (< 1 ,< 2 ,..,i B )2 4l « ia ...^, (4.76) 

n=l \/2/ i<j 1< j 2< ... < j n <7v 

where 

z . = e x-cit-x i0 = e x E~\ (i = 1, 2, TV). (4.7c) 

Furthermore, to compare the limiting form of u resulting from the peakon limit with the 
iV-peakon solution given by ref. 21), we shift the phase constant appropriately, so that 



nil a 2 



2^- (A, -A,) 2 A 2 ' 



(i = l,2,..,JV). (4.8) 



We substitute (4.8) into (4.7) and use (4.6) to modify them into the form 



/i~l + E en( ^ 1} (| J d ne nx DP_ n , (4.9a) 



71=1 



/ 2 ~ 1 + 2 £ e^) (^j n d n+1 e nx D^_ n , (4.96) 
where the positive coefficients d n are defined by 

d n = d n (t)= [f 1 ' (n = 1,2, ...,7V). (4.9c) 

Thus, equation corresponding to (3.5) becomes an algebraic equation of degree N + 1 

N+l 

£ e^" 1 ^/"-"* 1 + 0(e JV ( JV+1 ) +1 ) = 0. (4.10a) 

n=0 

Here, as in the 1-soliton case / = /2//1 and the coefficients /i n are defined by the relations 

h = 1, (4.106) 
/>„ = d^- 1 )* (fD%l n - 2D$_ n+1 ) , (n = 1,2,..., TV), (4.10c) 
h N+1 = - ni=1 Ai e . (4.10d) 
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To evaluate u using (3.7), we need the derivatives f t and f x . They are derived simply 
from (4.10) by differentiation. Explicitly 

Y^AT+l n ( n -l)7 fN-n+1 

ft v n=l , (4.11a) 

sr^N+l n (n-l)u fN-n+1 

fx m n=1 — • (4.116) 

En =0 e^-V(N-n + l)h n f"-" 

It follows from (4.2c) and (4.10) that 

h n , x = {n- l)h n + d n e nx D^_ n , (n = 1, 2, .., TV), (4.12a) 

h N+hx = Nh N+1 , (4.126) 

N 

K,t = -(£ C *K + d n e in ~ 1)x {fD^ - 2D<$_ n+1 , t ) , (n = 1, 2, N), (4.13a) 
i=i 

'ijv+i,* = - hN+1 ' (4.136) 

The expression of u in terms of / now follows from (3.7) and (4.10)-(4.13). After an 
elementary calculation, we find that 

(Ef=i c,)f N + Eli e^d n e^ (e*D<Sl n>t - 2D<$_ n +i, t ) f N ~ n 

u j- ^ m J - . (4.14) 

f N + Eli e^d n e^D^_ n f N - n 

The unknown / is obtained by solving the algebraic equation (4.10). The procedure for 
constructing solution can now be carried out straightforwardly. Here, we summarize the 
result. 

We seek the series solution of the form 

/ = / (0) +e 2 / ( 2 ) +e 4 / (4) + ^ ) (415) 

under the conditions f^> > (n = 0,2,4, ...). Note that the odd powers of e have been 
dropped in the above expansion. An inspection shows that this can be justified since 
eq. (4.10) includes only the even powers of e. See also the 1-peakon case in §3. We find 
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that the leading-order asymptotic of the positive solution of equation (4.10) is expressed 
simply as 

/ „ / (2n-2) e 2(n-l) = _^ = _ e 2(n-l) > (416fl) 

if the following inequalities hold for h n (n — 1,2, N + 1) 

hi > 0, h 2 > 0, hn-i > 0, /i„ < 0, /i n+ i < 0, h N+1 < 0. (4.166) 

To obtain (4.16a), we assume that / has a leading-order asymptotic of the form / ~ 
e 2mj(2m) and su b s itute this into (4.10). Then, (4.10) is expanded in powers of e as 

[h m {f {2rn) } N ~ m+1 + h m+1 {f {2m) } N ~ m ] e m(2JV " m+1) + (e m{2N - m+1)+2 ) = 0. (4.17) 

By taking the coefficient of e ™i™- m +i) zerQj we find that y(2m) = _ hm+i / hm ^ ). 
According to (4.16b) and the requirement /( 2m ) > 0, the integer m is determined uniquely 
as m = n — 1, which immediately leads to (4.16a). Note from (4.10c) that the coordinate 
x = x n giving rise to the equality h n = can be specified as 

l(0) 

(n = 1,2, ...,7V), (4.18) 



x„. = In 



D (2) 
U N~n 



and /i„ > for x > x n (n = 1,2, JV). An important consequence deduced from (4.18) 
is the notable inequalities 

Xj-i<Xj, (j = l,2,...,iV + l), (4.19) 

where we have used the convention xo = — oo and x^ + i = +oo. In fact, we substitute 
(4.18) into (4.19) to rewrite them in the following alternative forms 

D%l j+1 D%l j+1 - Df_ j+2 D%_ 3 > 0, (j = 1, 2, N + 1). (4.20) 

The left-hand side of (4.20) equals to (jD$_j + ^ by Jacobi's identity (4.5) with m = 
and n = N — j and consequently it always has a positive value, which proves (4.20). We 
see from (4.10c), (4.18) and (4.19) that under the inequalities (4.16b), x must lie in the 
interval < x < x n {n = 1, 2, N + 1). 



12 



Let u n be the waveform of u(x, t) for x in the interval x n _i < x < x n {n = 1, 2, N) 
at any instant t. We substitute (4.16) into (4.14) and see that both the denominator and 
numerator of (4.14) have a leading-order asymptotic of order e ( n - 1 )( 2N - n ) . Consequently 
the limit exists when e tends to zero. It turns out that u n has a limiting waveform given 
by 

u n = — -, (4.21a) 



where 



F n = -d n ^DP_ n+1 h n + d n e x D^_ n h n . x , (4.216) 

Gn = d n - ie -» (e x D^_ n +i,t ~ 2^ 0) _ n+2 ,t) h n - d n (fD%\ t - 2D^_ n+1 , t ) K-i- (4.21c) 
Inserting h n from (4.10c), F n becomes 

F n = 2d n . 1 d n e^ (4 n+ i4n + i - D { Z n D^_ n+2 ) . (4.22) 

We use the formula (4.5) with m = and n replaced by N — n in (4.22). Then, F n 
simplifies to 



F n — 2,d n —id n €^ ^ 
By a similar calculation, G n is transformed to 



(i) 

N-n+l 



(4.23) 



G n = 2d n . 1 d n e^ x 



e x r 

y i 



n (2) n (2) n (2) n (2) 

^N-n+lfi^N-n ^ N-n,^ N-n+l 



+D m D (o) , D (o) 



r,(2) _n( 2 ) nW _ nW ni 

N-n+l,^ N-n+l ^ N-n+lfi^ N-n+l 1J N-n+2,t ^ N -n 



l(0) 



l(0) 



i(2) 



+2e~ :C { 



D 



(0) 



(0) 



N-n+2,t ±y N-n+l 



-D 



(o) 



(0) 



N-n+l,t LJ N-n+2 



(4.24) 



The following identities among determinants are verified with the aid of the Jacobi identity 
(see Appendix B): 

>(2) 



u N-n+l,t u N-n 



n(°) n(°) /;"• // - , , 

^ N '-n+2,t^ N -71+ 1 - Ly 7V-n+l,t- Ly Af-n+2 — z - Ly 7V-n+2- Ly Af-ri+l ) 



(2) n (2) n (2) _ 9n (.ij n 



N-n+l^N-ni 



7-)(2) 7-)(0) _ D W 



(2) 

AT-n+l.t-^W-n+l 



N-n^N-n+l 

ziy A r -n+2 iy JV- 
,(0) nn(-l) /i.o\nW 



2-0^-1+2(1) 2)-Dat-?i+H 



(4.25a) 

(4.256) 
(4.25c) 
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n(°) D {2) - D (0) D {2) -2D { ~ 1) (]-2)D {1) 

^N-n+lJ^N-n+l ^ N-n+2,t LJ N-n ~ ZLJ N-n+2\ l i Z ') LJ N-n^V 

Substituting (4.25) into (4.24), we can reduce (4.24) considerably, giving rise to 



N-n+1 



e N—n + ^ e u N-n+2 



(4.26) 



The expression of u n now follows from (4.21a), (4.23) and (4.26). It is expressed compactly 
in terms of the Hankel determinants as 

u n = eDN ~ n+ { ^ Dn - + \ (n = l,2,...,iV). (4.27) 

D N-n+1 

Note that when n — 1, (4.27) reduces to u\ = e x D N ^_ l /D$ due to the relation = 

which represents the waveform of u in the range x < x 1 . Using (4.6), this expression can 
be written in the form 

N 

Ui = Y^ b ^ (4.28a) 
i=i 

with 

N N 

b= U A ? II ^ - X i)~ 2 K 3 , (i = 1, 2, .., JV), (4.286) 

j=l J'=l 

where Aj and Z; L are defined respectively by (4.2c) and (4.7c). 

For x in the range Xn < x, we find that w = Un+i has a particularly simple form 

N 

u N +i = 2e- x D?J = 2 CiZ~\ (4.29) 

i=l 

where we have used = A , (4.2b) and (4.7c). 
4-3 Comparison with the N-peakon solution 

Here, we compare the expression of the peakon limit arising from the iV-soliton solution 
with the iV-peakon solution given by Beals et al. 21 -* If we identify the quantities A™ (m > 
0) and A° introduced in ref. 21 with (m > 0) and Dn \ respectively, then we can 

write the latter solution in the form 

N 

u(x, t) = J2 )e~ |x ~ Xi(t)l , (4.30a) 
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rrii 



N-i+l 1 ^ N-i 
^N-i+l^N-i 



Xi = In 



2D {0) 

zlJ N-i+l 



D 



(2) 
N-i 



(t = l,2,...,N), 



(i = 1,2, ...,N). 



(4.306) 
(4.30c) 



When x lies in the interval x n -\ < x < x n , then u = u n takes the form 

n-l AT 



-(Xi-X) 



1=1 



JV-1 4 £) 



E 



»(0) \ 2 v-n fn(2) N2 

'i+i 



7V-n £) 



D (l) 
-^j+l-^i 



(4.31) 



The following formulas are particularly useful to simplify (4.31): 

2 



N—n ( JJ^ 

E " 



(3) 
AT— ra 



u i+\ u i 



D 



(0) 
i+1 



Af-1 

^(1) ^(1) 



(1) 

AT-n+1 

(-1) 
Af-n+2 



(i) 

AT-n+l 



(n = 0,1, ...,7V), 



(n = 2,3,..,7V + l). 



(4.32a) 



(4.326) 



To prove (4.32a), we write Jacobi's formula (4.5) with m = 1 and n = i — 1. It reads 



(4.33) 



Dividing (4.33) by D^D^ and summing up the resultant expression from % — 1 to 
i = N — n, we obtain (4.32a) by taking account of the relations = 1 (m = 1, 2, 3). 
By a similar calculation using (4.5) with m — — 1 and n — i, formula (4.32b) follows 
immediately upon noting the relation D^ N }\ = 0. If we introduce (4.32) into (4.31), we 
find that the resultant expression coincides with the peakon limit of the iV-soliton solution 
(4.27). 

Last, for xn < x, (4.30) becomes 



u = u N+1 = 4e x ^ 



N-l ( D W 



D (i) D {iy 



(4.34) 
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Using (4.32b) with n — N + 1, we can recast (4.34) to 

N 

u N+1 = Ae- X D { - 1] = 4e~ x ^ X' 1 ^, (4.35) 

i=i 

which is in agreement with (4.29) by (4.2c) and (4.7c). In conclusion, we have completed 
the proof that the iV-soliton solution of the CH equation reduces to the iV-peakon solution 
of the dispersionless CH equation in the peakon limit. 

5. Concluding Remarks 

We have developed a novel limiting procedure for recovering the nonanalytic iV-peakon 
solution from the analytic iV-soliton solution of the CH equation. In the process, the 
asymptotic analysis has been performed which expands all the wave parameters in powers 
of the small dispersion parameter. We have employed the key identity, i.e. Jacobi's 
formula (4.5) for determinants by which the significant part of the calculations has been 
carried out in quite a transparent manner. In view of this formula, we were also able to 
obtain the new formulas (4.27) and (4.28) which represent the waveform of the iV-peakon 
solution in a specified interval of the space variable. Our approach is quite different from 
that used by Beals et al 20 ' 21 ) which relies on the classical moment problem originally 
studied by Stieltjes. In a future work, we have a plan to apply our method to the N- 
soliton solution 27,28 ) of the Degasperis-Procesi (DP) equation to show that it recovers 
the iV-peakon solution already given by Lundmark and Szmigielski. 29,30 -* While both the 
CH and DP equations have a different mathematical structure, the iV-soliton solution of 
the DP equation can be written in a parametric form just like (2.1) and (2.2), the only 
difference being the structure of the tau-functions f\ and ji- This observation will make 
it possible to perform the peakon limit of the iV-soliton solution in a manner similar to 
that developed here for the CH equation. 
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Appendix A: Proof of (4.6) 

We introduce the n x N matrix P and the N x n matrix Q 



P = 



Q 



/ 1 1 • 


• 1 \ 






Ai A2 


AjV 




(A 


1 


\n— 1 \n— 1 


N ) 






?i X[ n+1 E 1 ■ ■ 


ym+n- 


' 1 E l \ 




71 \m+l 771 
1/2 A 2 -C/2 


ym+n- 




(A 









\\%E N X™ +1 E 



N 



\T n - 1 E N ) 



We calculate det(PQ) by two ways. First, multiplying P and Q and using (4.2), we obtain 

det(PQ) = £><">. (A -3) 

On the other hand, the determinant is expanded by means of the Binet-Cauchy formula 



as 



31) 





1 


1 


1 


det(PQ) = J2 




A« 2 


A in 


l<i 1 <i2<-<i n <N 


\ 71—1 


\ n—l 

K 





\ m 771 \m-t-± 77 
\mrp \n+lp 



\m 77 \m+l 77 



^ A„(ii,i 2 , ...,i n )(Xi 1 Xi 2 ...Xi n ) m Ei 1 E i2 ...E in . 

1<U<1 2 <. ..<Jn<iV 



ym+n-l 

K +n - l E- 2 



\f +n - l E,, 

(A-4) 



In passing to the last line of (A-4), we have extracted the factor \™E i{ from the lih row 
(I = 1,2, ...,n) and used (4.1). Formula (4.6) follows immediately from (A-3) and (A-4). 

Appendix B: Formulas for Determinants 

Here, we prove (4.25). We define the quantity D by the relation 



n _ n( 2 ) n(2) _ r>(2) r>(2) 



(P-l) 



We differentiate d! 2 -* by t with use of the relation dA m /dt = 2A m _i which is a consequence 
of (4.2b). The differentiation of the jth column of Dn^ is then proportional to the (j — l)th 
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column. Thus, it vanishes identically except for j — 1. Hence 



D% = 2 



A 3 ••• A n+1 
A2 A 4 ■ ■ ■ A n+ 2 



(B-2) 



A n A n+ 2 ' ' ' A2 n 

Using (B-2) and the relations (4.4), (B-l) becomes 
£> = 2^ 2 (n + 2;2)I^ 2 (n + l,n + 2;^^ 

It follows by using Jacobi's formula (4.3) and a fact -D^+ 2 > ( see (4-6)) that 

4 1 ] 2 (n + l,n + 2;l,n + 2) 



B-3) 



D 



(i) 

n+2 



£>2 2 (n + 1; l)I># 2 (n + 2; n + 2) - Z>£ 2 (n + 1; n + 2)D£ 2 (n + 2; 1) 



)(!) 



(1) 



>(!) 



(5-4) 



D2 2 (n+l,n + 2;2,n + 2) 



(i) 

n+2 



£>£J 2 (n + 1; 2) J D„ i / 2 (n + 2; n + 2) - D™ 2 (n + 1; n + 2)2?£ 2 (n + 2; 2) 



)(!) 



(1) 



(S-5) 



Substituting (B-4) and (B-5) into (B-3), we obtain 
2D%(n + 2;n + 2) 



D 



D 



(i) 

n+2 



£>£> 2 (n + 1; l)D%(n + 1; 2) - £>£ 2 (n + 1; 2) J D« 2 (n + 1; 1 

(5-6) 

We apply Jacobi's formula to the right-hand side of (B-6) and use the relations (4.4). 
Then, (B-6) is recast to a simple form 



D = 2D^ +l D^. 



(B-7) 



Formula (4.24a) is a consequence of (B-l) and (B-7) with n replaced by iV — n. Formula 
(4.24b) follows from (4.24a). To show this, we may simply replace n by n — 1 and Ei by 
X~ 2 Ei(i = 1,2, ...,7V), respectively in (4.24a). 
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By repeating the similar calculation, we can show that 

D®D<® 2 - D^D^l = -2d£>(1, n + 3; 2, n + 3)£>£S(n + 2, n + 3; 1, 2) 

= -2£>tS(l;2)I>£ 1 . (B-8) 

If we replace n by N — n in (B-8), we obtain (4.25c). Formula (4.25d) is derived by the 
same way. 
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